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Abstract 

What is the simplest Hamiltonian which can implement quantum 
computation without requiring any control operations during the com- 
putation process? In a previous paper we have constructed a 10-local 
finite-range interaction among qubits on a 2D lattice having this prop- 
erty. Here we show that pair-interactions among qutrits on a 2D lattice 
are sufficient, too, and can also implement an ergodic computer where 
the result can be read out from the time average state after some 
post-selection with high success probability. 

Two of the 3 qutrit states are given by the two levels of a spin-1/2 
particle located at a specific lattice site, the third state is its absence. 
Usual hopping terms together with an attractive force among adja- 
cent particles induce a coupled quantum walk where the particle spins 
are subjected to spatially inhomogeneous interactions implementing 
holonomic quantum computing. The holonomic method ensures that 
the implemented circuit does not depend on the time needed for the 
walk. 

Even though the implementation of the required type of spin-spin 
interactions is currently unclear, the model shows that quite simple 
Hamiltonians are powerful enough to allow for universal quantum com- 
puting in a closed physical system. 
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1 Introduction 



To understand which sets of quantum control operations are sufficient for 
quantum computing is still an important issue of research. Whereas the 
standard model of the quantum computer is based on one- and two-qubit 
unitaries there is meanwhile a large number of alternative proposals, e.g. 
computing by measurements only [TJ [2] or adiabatic computing [3l [U [5], El [7] . 
The latter model encodes a computational problem into an interaction such 
that the ground state of the Hamiltonian indicates the solution. In order 
to understand the power of adiabatic computing it is particularly interest- 
ing to know to what extent the transition into the ground state of simple 
interactions can already be a non-trivial computation process. Since it has 
been shown that every problem in the complexity class QMA can be encoded 
in a nearest-neighbor interaction of qubits located on a 2D lattice [8] it is 
clear that, indeed, relatively simple Hamiltonians are sufficient for adiabatic 
computing. This highlights the computational power of simple Hamiltoni- 
ans from one point of view. Another aspect of the "computational power of 
Hamiltonians" has been studied in models where the natural time evolution 
is interspersed by external control operations (see e.g. [TUl [TT], [12]). The 
intention of the present paper is to understand to what extent models with 
simple interactions can already be autonomous devices with full quantum 
computing power in the sense that quantum computation reduces to the fol- 
lowing protocol: (1) prepare an initial state in the computational basis which 
contains the data, (2) wait for a sufficiently long time, and (3) measure a 
sufficiently large set of qubits in the computational basis. In [13] we have 
constructed a Hamiltonian satisfying these conditions with the additional 
feature that the readout need not necessarily be performed within a specific 
time interval; in our "ergodic quantum computer" the result is also present in 
the time average of the dynamics after some simple post-selection with high 
success probability. Even though this property may be of minor practical rel- 
evance, we considered it as a necessary feature of an autonomous computer 
since one would otherwise require a clock as an additional device. Here we 
construct an interaction that is even simpler than the one constructed in [T3] 
since the latter requires 10-qubit interactions. 

The motivation to consider autonomous quantum computers is twofold. 
First, it could trigger new ideas how to reduce the set of necessary control 
operations in current implementation proposals by using the "natural power" 
of the interactions. In addressing this issue, the ergodic model defines only 
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an extreme case; realistic perspectives for quantum computing could arise by 
combining it with more conventional approaches. Second it is an interest- 
ing fundamental issue in the thermodynamics of computation how to realize 
computation in closed physical systems. Benioff, Feynman, and Margolus 
have already presented such Hamiltonian computers [T^ , [T5l [16] . Margolus' 
asynchronous Hamiltonian cellular automaton (CA) has the appealing fea- 
ture to be lattice-symmetric. In other words, it is driven by a finite-range 
interaction among qubits3 located on a 2D lattice such that the total Hamil- 
tonian is invariant under discrete translations. The synchronization of the 
CA is realized by a kind of spin wave propagating along the lattice and trig- 
gering the update of the cells according to some computationally universal 
update rules which are not specified any further in [TO]. However, its clock 
wave has to start in an uncertain position in order to obtain a well localized 
momentum distribution with mainly positive momenta. A localized wave 
front would also propagate backwards and would therefore not trigger a cor- 
rect computation. In [13] we have chosen the same synchronization scheme 
but we start with a localized wave front since we allow only for preparations 
of basis states. The fact that the dispersion of the wave front leads to com- 
pletely undefined computation steps is irrelevant since the time average of 
the dynamics encodes the correct result. The feature of our ergodic model to 
show the correct computation result also in the time average was hence only 
a nice byproduct of the fact that we must use a concept which works with 
a strange kind of clock: The latter starts with a well-defined time but then 
it counts backwards and forward with completely undefined counting speed. 
Note that the time average can also be considered as the "generalized final 
state" of our computer since a final output state in the usual sense cannot 
exist for finite dimensional Hamiltonian models. 

To ask for a simple finite range Hamiltonian that is universal for quantum 
computation is in some sense similar to asking for a simple computationally 
universal quantum cellular automaton (for some proposals see [HI [19]) with 
the only difference of considering update rules which change the state only 
in an "infinitesimal" way. However, local interactions in lattices have typi- 
cally the property to spread the information into increasingly large regions, 
whereas it is possible to construct update rules for cellular automata that 
work by propagating the information forward column by column [18] . This 

1 In [T7] it is furthermore argued that one of this clock Hamiltonians for CAs in one 
dimension is quite close to real solid state interactions. 



3 



apparent difference between discrete and continuous time dynamics can al- 
ready be explained with the translation operator on n qubits: the cyclic shift 
is a unitary that can be achieved by local update rules [20], whereas the 
Hamiltonian obtained from the logarithm of the shift contains interactions 
between distant qubits. Hamiltonians with finite interaction length will typ- 
ically spread the information over the lattice. Therefore, one has to be more 
modest and demand that the correct result can only be present with high 
probability. 

The paper is organized as follows. In section [2] we explain the hardware 
of our model consisting of a chain of atoms on a 2D square lattice. We 
construct a nearest-neighbor interaction which leads via a simple effective 
Hamiltonian to a coupled random walk that will later trigger the implemen- 
tation of gates acting on the atom spins. In section [3] we explain how to 
arrange spin-spin interactions among the atoms such that they implement 
logical gates based on holonomic quantum computing. In section H] we de- 
scribe the complete Hamiltonian of our computing device. In section [5] we 
show that the time evolution of the effective Hamiltonian is solvable since it 
can be transformed into a quasifree evolution of fermions. Based on this solu- 
tion, we estimate the time required for the computation process. In section [6] 
we briefly sketch some ideas on the realization. Section [7] describes how to 
construct the Hamiltonian such that it implements a universal quantum cel- 
lular automaton. This is to obtain programmable hardware. In the appendix 
we present an alternative model for the atom propagation which is caused 
by an even simpler Hamiltonian since the nearest-neighbor interactions in 
the model presented in the main part involves also diagonal neighbors in the 
lattice. In contrast, the model in the appendix uses only nearest neighbors in 
a strict sense. However, the corresponding Hamiltonian dynamics will not be 
solvable. We will therefore consider an incoherent analogue and show that 
the corresponding classical random walk would trigger the implementation 
of gates in the desired way. For the coherent model, we can only conjecture 
that an appropriate propagation is achieved. 

2 The synchronization Hamiltonian 

Now we present a model of a 2D lattice with a coupled propagation of a 
chain of atoms which will later be the synchronization mechanism of the 
computer. While propagating along the horizontal direction (see Fig. [T]) the 
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Figure 1: Initial atom configuration for a 7 x 7 lattice. The digits at the left 
indicate the corresponding binary string which will be introduced in section [5] and 
the arrow shows the only possible motion. 



spin of the atoms (which represent the logical states) will be subjected to 
spatially inhomogeneous interactions that implement the desired gates. It is 
important that the chain does not tear since the spin-spin interactions to be 
described later will only be active between adjacent atoms. Furthermore, it 
is important for the synchronization that the chain remains connected. This 
will become clear in section [3J 

We will now describe the synchronization Hamiltonian in detail and show 
that it generates a diagonalizable quantum walk. The lattice and the initial 
atom configuration are shown in Fig. [IJ The lattice sites are given bynxn 
black fields of the modified chess-board in Fig. [IJ The 2n — 1 atoms are 
only allowed to move forward or backward along the horizontal direction, 
i.e., along the rows. Each atom can only hop from a black site to either 
the next or the previous black site in the row. Furthermore, the atoms can 
only move such that the chain does not tear, i.e., atoms in adjacent rows are 
always diagonal neighbors. Fig. [2] shows a possible configuration. 

Now we describe the Hamiltonian that allows only collective motions of 
the atoms respecting these rules. The grid of black fields is labeled by 
with i,j = 1, . . . , n. The sites with j = 1, . . . , n indicate the diagonal 

line leading from the leftmost field to the uppermost one and sites (i, 1) with 
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Figure 2: Possible atom configuration with its corresponding binary string (as 
explained in section ED and the possible motions. 



i = 1, ... ,n indicate the diagonal line from the leftmost field to the lower- 
most one. We consider the Hilbert space of all atom configurations as a 
subspace of an n 2 -qubit register of the form 

n q : = (c 2 )®" 2 , 

where qubit corresponds to lattice site and the basis states |0) and 
|1) indicate that there is no atom or that there is an atom, respectively, at 
position A subspace of 7i q is the clock Hilbert space H c spanned by 

all allowed configurations. They correspond to those basis states in 7i q for 
which there is exactly one atom in each row and for which the atom chain is 
connected. First we introduce interactions between site (i, j) and (z + 1, j + 
which generates independent hopping of atoms along the rows: 

n-l 

K := Yl ( a M' a l+i,j+i + h.c^j . (1) 

To achieve that the atom configuration remains a connected chain we intro- 
duce a strong attractive force between diagonal neighbors: 

H pot :=-E NijN k>l + E l, (2) 

<(i,j),(k,l)> 
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where E G M + is assumed to be much larger than 1 and Nij = aJ.-Oij is the 
projection on the states with an atom in position The sum runs over all 
unordered pairs of sites (k, I) with \i — k\ + \ j — l\ — 1. The physically 
irrelevant term Eq is (for purely technical reasons) chosen such that the 
initial configuration has zero energy. The energy is minimal if all attractive 
interactions are active. Otherwise, if the chain is not connected, the potential 
energy is at least E. We define the "synchronization Hamiltonian" , which 
leads to the desired coupled motion of atoms, by 

H:=K + H pot . (3) 

In order to analyze the dynamics generated by H we first argue that it can be 
replaced with an effective clock Hamiltonian H e ff having 7i c as an invariant 
subspace. The initial atom configuration defines a ground state of H pot . Due 
to E ^> 1 the comparably small perturbation K cannot circumvent the gap 
to one of the first excited states of H pot . If P denotes the spectral projection 
of Hp t corresponding to the ground state energy we will therefore expect 
a dynamical evolution according to 

H eff :=PKP. (4) 

Note that not all states in the image of P are allowed clock states since 
the image of P contains also states with more than one atom in a row. 
However, these ground states could only be reached from an allowed state if 
K contained hopping terms along columns. 

The following lemma (proven in the appendix) justifies more formally 
that we may analyze the dynamics according to H e ff instead of H: 

Lemma 1 For E > n 6 and n > 10 the norm distance between the effective 
Hamiltonian in equation §Jj§ and the restriction He /2 of H to the eigenspace 
corresponding to eigenvalues smaller than E satisfies 

n 3 

\\He/2 — H e ff\\ — ■ 

If we increase the attractive part of the interaction proportional to n 6 
if n increases the time evolution induced by H e ff is a good approximation 
for the true one. If the intention of this article was to propose a physical 
implementation scheme this would be a rather bad scaling, leading possibly 
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to hard practical problems. However, from the computer science point of 
view, it is essential that the increase of energy is only polynomial in n and 
we will thus obtain an efficient model of computation. 

We will now give a more explicit form of H e ff. To be more precise, we 
will describe an operator whose restriction to 7i c coincides with H e ff. Let |c) 
denote some basis state in TC C , i.e., c describes an allowed atom configuration. 
Then 

n-1 

H eff \c) = PKP\c) =PK\c)=PY, ("<,",.,,.: + h.c) \c) (5) 

is a superposition of all atom configurations with connected chain that can 
be reached from c by moving one atom forward or backwards in its diagonal 
row. In other words, each hopping term o^aj+j is only active if the sites 
+ l and i + 1, j are occupied. We have therefore a conditional hopping 
given by 

n-1 

H eff = Yl a id a t+l,j+l N i,j+l N i+l,3 + h - C - w > 

where the symbol TC C at the right indicates the restriction to TC C . Using a 
slightly more intuitive notation, which respects the orientation of Fig. [TJ the 
summands of H e ff can be denoted by 

N 

a + h.c. , (6) 

where this 4-local operator acts on 4 black sites enclosing a common white 
site. 

We are interested in the dynamical evolution of the quantum state defined 
by the initial atom configuration. Since it is a ground state of the dominating 
Hamiltonian H the major part of the state vector is contained in that spectral 
subspace of H which corresponds to energy values not greater than E/2. 
Using the bound of Lemma[T]we may then estimate the error in the dynamical 
evolution caused by replacing the true evolution with the evolution according 

tO H eff . 
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Theorem 1 Let G 7i c be an arbitrary allowed atom configuration. Then 
we have 

||( e -^_ e -^//*)|^)|| <et + 2n ' 2 

where 



n 

e := 9 



. E 

3 



E 



as in LemmaUl 



Proof: Due to H pot \ip) = the average energy of \ip) satisfies 

(tp\H\ip) = (if)\K\ij;) < \\K\\ < n 2 , 

where the last inequality is given by counting the number of terms in eq. (PQ). 
Let Q be the spectral projection of H corresponding to all eigenvalues less 
than E. Then we have 

(ip\QHQ\ip} + - Q)H(1 - QM) = (ij>\H\1>) < n 2 . 

Using 

^|(i-g)i/(i-g)|^)> J E;||(i-g)^)|| 2 , 

and 

(iP\QHQ\iP) = {ii)\K\if)) > -\\K\\ > -n 2 

we conclude 

||(l-Q)|^|| 2 <^. (7) 

We have 

|| (g-ifft _ giffe//*) |^|| < ||( e -^/ 2 *_ e -^//*)Q|^)|| 

+\\(e- i&t -e- iH «ff t ){l-Q)\il>)\\ 
< et + 2\\{l-Qm\\, (8) 

where we have used || exp(iA) — exp(iB)\\ < \\A — B\ for two self-adjoint 
matrices A, B. The statement follows then by replacing the last term in 
eq. (jHJ) with the right hand side of ineq. ((7j). □ 



One could use the condition t <C vE/(9n 3 ) of Theorem [T] to determine 
the time scale for which H e ff is a good approximation when E is given. We 



9 



can also assume that the time scale is given and we have to determine the 
interaction strength E. After all, we have to ensure that the approximation 
is valid for all t < to where to is some a priori given upper bound on the time 
required by the clock wave to pass the relevant region on the lattice, i.e., the 
region containing the spin-spin interactions which will be described in the 
next section. 

The effective Hamiltonian leads to a quantum walk on the space of atom 
configurations which will be diagonalized in section [51 In the next section we 
will explain how this walk may trigger the implementation of gates. 

3 Holonomic implementation of logical gates 

Now we consider not only the positions of the atoms but also their inner 
degree of freedom, e.g. their spin. We replace the qubit at each lattice site 
with a qutrit and extend hence the space 7i q to 

n := (c 3 ) 0n2 . 

The basis states |0), | j), I T) indicate the absence of the atom or its two 
possible spin states, respectively. Before we describe how to add a spin-spin 
interaction Hamiltonian to H which leads to the implementation of gates 
when the atoms are moving along the rows, we should first mention two 
obvious approaches to imprint gates by interactions and why they are not 
suitable for our goal. 

In [13] and (16] the gates are directly coupled (as additional tensor product 
components) to the synchronization Hamiltonian. Certainly we could define 
the hopping term in eq. (pQ) between column j and j + 1 in any desired 
row such that the spin is not conserved during the atom propagation but 
subjected to some rotation. This would provide us with one-qubit gates. 
Such a one-qubit gate would clearly be inverted when the atom moves back 
to column j again. The conditional state of the considered qubit, given that 
the atom is found on the right of column j would hence be subjected to the 
desired transformation, no matter whether the atom has traveled back and 
forth several times as is is typical for a random walk. 

To imprint two-qubit gates into the Hamiltonian is more difficult. The 
obvious method to couple their implementation with some coupled motion 
of two adjacent atoms would require more non-local interactions than only 
two-particle terms. The second obvious method would be to complete the 
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clock Hamiltonian by an additive spin-spin interaction term between some 
adjacent sites in the same column. This leads to the following problems. 
First the atoms do not stay there for a well-defined time. Some part of the 
wave packet moves already and one part stays. Second the atoms travel back 
and forth and pass the interaction region several times. Both effects would 
in general entangle atom position and logical spin states in an uncontrollable 
way. 

A solution to the latter problem is holonomic quantum computing [2TII22] . 
The idea of this approach, in the usual setting, is that some time-dependent 
Hamiltonian G(t) is adiabatically changed along a closed loop such that, 
on an appropriate degenerate subspace, the overall effect is a unitary which 
depends only on the loop (described by representing G(t) in some parameter 
space) and not on the speed of the change of G(t). The unitary in the end 
remains the same even if the vector G(t) has been moving back and forth on 
this loop in the parameter space. 

We apply this concept to our model where we achieve the time-dependence 
of the spin-spin interaction by the motion of the atoms when they pass spa- 
tially inhomogeneous interactions. We imprint a spin-spin interaction (vary- 
ing slowly along the rows) such that is describes a closed loop after the atoms 
have passed a certain region. Then the spin state of the atoms does not de- 
pend on the time required for the passing, neither does it depend on the 
number of times the atoms were traveling back and forth before the whole 
region has been passed. This ensures that the entanglement between the 
position degree of freedom and the spin is lost as soon as the region has been 
passed by the whole atom chain. Before we explain how to make use of this 
idea in detail, we rephrase a result in section II of Ref. [23], which is useful 
for US, clS db lemma: 

Lemma 2 (Holonomic Gates) 

Given a family of Hamiltonians G(t) on a finite dimensional Hilbert space 



where X is a self-adjoint operator on Ji such that the family G(t) is a closed 
loop, i.e., G(0) = G{T). Let the change of G(t) be sufficiently slowly to 
consider it as an adiabatic change, i.e., for each entry Gij(t) of G(t) we 
have 



H by 



G(t) = e tXt G e 



-iXt 



with t e [0, T] , 



d 
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where A is the smallest spectral gap of G(0). 

Let Ut be the time evolution generated by the family G(t) after the closed 
loop. Let /C be an eigenspace of Go with eigenvalue X. Then Ut keeps JC 
invariant and its restriction is given by 

e - t XT e iQXQT^ (g) 

where Q is the projection onto /C. The unitary (TJ|) is also called a non-abelian 
geometric phase. 

Since the holonomic concept implements the well defined unitary 

exp(iQXQT) 

only on a subspace K, of the full Hilbert space our 2n — 1 spin particles cannot 
be used for 2n — 1 logical qubits. Instead, the atom in row 2i — 1 will encode 
a logical qubit together with the atom in row 21 Moreover, we can only 
use rows being not too far from the row in the middle since the others are 
too short to imprint interactions that implement gates. Two adjacent spin 
particles encode one qubit with logical states \0)i, \ via 

\0)l := I I) ® I T) and |1>, := | |) ® | |> . 
We will call their span C the "code space" . 

One-qubit gates 

To explain how to implement one-qubit gates we first restrict our attention to 
two adjacent rows. We add interactions as follows. The whole region which 
carries the gate consists of a stripe of width 2 as shown in Fig. El We call this 
region "interaction stripe" . Inside this stripe, we add a "gate Hamiltonian" 
to H which consists of spin-spin interactions Vj between certain diagonal 
neighbors in adjacent rows (as indicated by the edges labeled with V\, . . . , VJ 
in Fig. [3]). Each Vj is a pair-interaction between two sites. We choose the 
following interactions: 

■ 27r(j— 11 v -2^(1-1) v 

V j :=e l ^ ±x (o z (^N + N®a z )e- l ^ lx , (10) 

with j = 1, . . . , I. Here X is some self-adjoint operator on C 3 ® C 3 specified 
later and 

W:=U)UI + lTXTl 
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Figure 3: Spin-spin interactions for a logical one-qubit gate. Note that these are 
nearest neighbor interactions in a 2D lattice when the grid is oriented in diagonal 
direction with respect to the pattern of the chess-board. 



is the particle number operator in analogy to the operator Nij in eq. (j2J). 
The Pauli matrix a z has to be read as formally acting on the space C 3 of 
the corresponding qutrit, even though it is zero for the state |0) and acts 
therefore only on the spin states. It is important to note that the degenerate 
eigenspace of Vo and V/_i coincides with the code space C. If the spins are 
initialized such that their state is in C before they pass the interaction region 
they will at every time instant remain in the degenerate subspace of the 
interaction that is active (provided that the change between different Vj is 
adiabatic). This will implement the "non-abelian phase" exp(iQXQT) of 
Lemma EJ 

The length I of the interaction stripe is chosen such that the change of 
interactions from j to j + 1 can be considered as approximating a continuous 
change and furthermore as adiabatic when the atom chain is propagating in 
horizontal direction. To sketch how to choose I we define a typical time scale 
t on which an atom jumps from one site to its neighbor. Then we have to 
ensure that each entry V- k of Vj satisfies 



V j+1 V j 



3 



«A, (11) 



where A is the smallest energy gap of V\. Here, the time r is defined as a 
dimensionless quantity of the order 1 since the hopping terms aija i+1 j +1 +h.c. 
in eq. (CQ) appears with constant one. Physical dimensions are irrelevant since 
inequality (TTTT) is invariant with respect to a simultaneous rescaling of all Vj 
and the hopping terms in eq. (OQ). We obtain therefore the dimensionless 
condition I ^> 1. 
Now we may set 

X := o x (g> (cos (f)<J x + sin <j)<J z ) , (12) 
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with arbitrary angle 0, or 

X := a y (g> (cos (pa x + sin 0cr 2 ) , (13) 

where we have again slightly abused notation since we did not explicitly 
indicate the embedding of the 2-qubit operators in eqs. ffl2|) and ffl3|) (acting 
on the 4 possible spin configurations) into the two qutrit space. The following 
lemma shows which gates are implemented by the above choices for X. 

Lemma 3 Let two atoms pass a region where they are subjected to the Hamil- 
tonians V\, . . . ,Vi in eq. /[Tty) with X as in eq. / fjgj) or eq. ( T73jj . Assume 
furthermore that I is sufficiently large to consider the change of interactions 
as adiabatic. Then their encoded qubit is, up to an irrelevant global phase, 
subjected to the transformation exp(27r? cos <f)cr x ) or exp(2ir i cos <f) o~ y ) , re- 
spectively. 

Proof: On the space C 2 ® C 2 spanned by the 4 spin states the interaction 
changes according to 

■ 2-rrlj — 1) .2ir(j — 1) v 

Vj := e"^=r lx (a z ® 1 + 1 <g> a z ) e - l ^^ x , 

with j = 1, . . . , I. Due to Lemma [2] the code space C is then subjected 
to exp(i27iA) with A := PcXPc, where Pq is the projection onto C. For 
both choices for X, the relevant operator A consists then only of transitions 
between | f j) and | It), which corresponds to a a x if the transition amplitudes 
are 1 (as it is the case given by eq. (1T21 ) and to a y if they are i and —i (as 
in the case of eq. [13]) . □ 

Lemma [3] shows that we can generate arbitrary one-qubit transforma- 
tions by concatenating interactions with X as in eq. (fT2j) and (jTBl using 
appropriate angles 0. 

Two-qubit gates 

For the implementation of logical two-qubit gates we consider interaction 
stripes that consist of 3 adjacent rows (see Fig. H]), where row 1 and 2 belong 
to the first logical qubit and row 3 is part of the second logical qubit. We 
define interactions V and Uj such that V connects row 1 and 2 and Uj 
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Figure 4: Spin-spin interactions for a logical two-qubit gate. Row 1 and row 2 
encode one logical qubit and row 3 is part of a second qubit (together with row 4 
which is not shown in the figure). 

connects row 2 and 3 as shown in Fig. HI We will refer to the atoms in row 
1, 2, 3 as atoms 1, 2, 3, respectively. 

The interaction between row 1 and the row 2 is always V in the whole 
interaction stripe. Row 2 interacts with row 3 according to U%, U2 ■ ■ ■ , U\. 

As soon as atom 2 has entered the interaction stripe (shown in Fig. [4] for 
I — 7) it is simultaneously subjected to interactions V and U\. Moreover, V 
and Ui are turned off simultaneously as soon as atom 2 leaves the interaction 
stripe. The total interaction that is active on the system consisting of the 
spins of atoms 1,2,3 is therefore either 

V®l + l®Uj with j = 1, . . .,1, 

when atom 2 is inside the interaction stripe, or otherwise. 
We choose 

V:=a z ®N, (14) 

where the first tensor component refers to row 1 in Fig. [4] and the second to 
row 2. The operator N has no effect on the spin of the atom in row 2. It 
only ensures that the ^-Hamiltonian is not switched on before the atom in 
row 2 has entered the left-most black field in row 3 in Fig. [4} 
The interaction Uj changes according to 

U 3 := o z ® I I) (i I + (e'^^e- 1 ^*) ® | !><T | , (15) 

where X is an operator that acts on the one-qutrit space. It will be specified 
later. On the space C 2 <E> C 2 ® C 2 spanned by all 8 possible spin states, we 
have 

V®l+l®Uj = ff 2 ®l®l+l®a 2 ®| |) (| \+l®(e l2 - 1 ^ llx a z e~ i ^ llx )®\ |)(T I 
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The idea is that the adiabatic change of the degenerate Hamiltonian on 
row 1 and 2 is controlled by the spin of the atom in row 3, i.e., by the logical 
state of the second qubit. Whenever the state in row 3 is | f), rows 1 and 2 
are subjected to the Hamiltonians 

, 2^0-1) y ' 2Tr(j-l) y 

a z ®l + l®{e l «-i x a z e «-i x ) , 

otherwise they are subjected to a constant Hamiltonian. We find: 

Lemma 4 Let V and Uj be as in eq. [Lfy and eq. l[Po\) . respectively, such 
that they connect row 1 and 2 as well as 2 and 3, respectively. Choose the 
corresponding operator X as 

X := cos<p o~ x + sin o z , 

with arbitrary angle <fi. Consider the logical 2-qubit space C ® C, where the 
first tensor component refers to the atoms in rows 1 and 2 and the second 
to those in 3 and the additional row 4 (which is not depicted in Fig. [^). 
Then the adiabatic change of the interaction V (g> 1 + 1 ® Uj from j = 1 to 
j = I implements a controlled-exp(i sin <fr a z ) transformation on C <8> C with 
the lower qubit as control qubit and the upper as target (in the tensor product 
in eq. (F73J) this correponds to the right and the left side, respectively) . 

Proof: Let the second qubit be in the logical state |0)^. This means that 
the spins of atom 3 and 4 are in the state | jf)- However, for the interaction 
only the state of atom 3 matters. Since it is | |), atoms 1 and 2 are constantly 
subjected to the spin Hamiltonian 

a z <g> 1 + 1 ® a z . (16) 

Since the kernel of this operator coincides with the code space the logical 
state is not affected at all. 

Assume now that the second qubit is in the state |1) and atom 3 is there- 
fore in the state | |). Then atom 1 and 2 are subjected to the Hamiltonians 

, - 2^0-1) y . 27TQ-1) y 

a z ® 1 + 1® (e l '-1 x a z e 1 «-i x ) , 

with j = 1, . . . , I. Due to Lemma [2] the unitary that is implemented on the 
code space after such an adiabatic change is given by exp(i2irA) with 

A := P c (l ® X)P C = sin0P c (l <g> a z )P c . 
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Describing this operator with respect to the logical basis |0)z, \ we obtain 
the operator sin0<j z . This shows that the unitary exp(i2ir sm<j)<j z ) is imple- 
mented after the atoms have passed the interaction stripe. Since we have also 
shown that this unitary is only implemented if the state of the second qubit is 
we have know shown that the net effect is a controlled-exp(«27r sin (pcr z ) 
gate. □. 

We can imprint interactions like those above between any adjacent log- 
ical qubits as well as we may obtain one-qubit rotations as in Lemma [3] by 
appropriate interactions. It is known that conditional phase gates together 
with the set of one-qubit rotations allow for universal quantum computing 
|24j . Hence we are able to design interactions for arbitrary quantum circuits. 

4 The complete Hamiltonian 

We obtain the complete Hamiltonian H of our autonomous computing device 
as a sum of the synchronization Hamiltonian and the spatially inhomogeneous 
spin-spin interaction that implement holonomic transformations. It reads: 

H:=H pot + K+ W w',f > 

<i,j,i',j'> 

where H pot is the strong attractive force defined in eq. [2] and K describes the 
hopping terms 

71-1 
M=l 

where we have, in slight abuse of notation, adapted the definition of eq. for 
spin-less creation and annihilation operators into those for spin 1/2 particles. 
For instance, a\- j creates a particle with spin down. Similarly, the term H pot 
is adapted to spin 1/2 particles in the sense that the particle number operator 
Nij is given by 

The spin-spin interactions Wijyj> are only non-zero if site is adjacent 
to site (i',f) and if the pair of sites belongs to some interaction stripe. 
Inside these stripes, they are given by the interactions V,Vj,Uj described 
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L_ J output region 

IJ 

Figure 5: All spin-spin interactions implementing some circuit U are confined to 
the "circuit region" R. As soon as all atoms have left this region, the implemen- 
tation of U is complete. 

in the preceding section. We assume that all these rectangles ("interaction 
stripes") are confined to a square of length k, i.e., the spin-spin interaction is 
only non-zero for < k for some k < n. We call this region the circuit 

region (see Fig. [5]). As soon as all atoms have passed it, the whole circuit is 
implemented. The complement of the interaction region will be called output 
region since we will read out the result of the computation there. 

5 Diagonalization of the quantum walk 

Now we show that the Hamiltonian H leads to a propagation of atoms that 
implements the gates in correct order. First of all, we observe that the 
gates are irrelevant for the quantum walk. Different particle configurations 
correspond to mutually orthogonal vectors in the Hilbert space. Whether 
or not the particles are subjected to an additional change of their inner 
degree of freedom is irrelevant for the dynamics as long as we consider the 
adiabatic limit where the effect of the spin-spin interaction is only that it 
implements unitary gates on the spin states. Hence we will analyze the 
dynamics generated by H e ff instead of the dynamics generated by H to derive 
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the relevant probabilities for finding the atom chain at certain positions. 
After recalling that the adiabatic change of the spin-spin interactions has 
implemented the desired gates given that the atom chain has passed the 
interaction stripes we have then shown that the circuit has been implemented. 

To diagonalize the quantum walk we show that it is isomorphic to a 
walk of free fermions propagating on a one-dimensional chain. On the left 
of Figs. [T] and [2] we have characterized the configurations by binary words of 
length 2n — 2 =: 2m. The symbols 0, 1 as jth digit indicate whether the atom 
in row j + 1 is behind the atom on row j or in front of it, respectively. The 
initial configuration is hence characterized by m = n — 1 symbols followed 
by m symbols 1. Since the atom on the very top and that one on the very 
bottom are fixed, the number of symbols 1 remains constant and the vector 
space spanned by the possible atom configurations is therefore the subspace 
Ti, m of (C 2 )® 2m spanned by words with Hamming weight m. With respect to 
such a representation, the effective clock Hamiltonian H e ff, when restricted 
to TC C , consists of operators which replace some pattern 10 by 01 or vice 
versa. It can be written as 

2m- 1 2m- 1 



H s := l^ 1 ®|0)(l|®|l)(0|®l $52m -- ? '- 1 + /i.c. = ^2b}b j+1 + h.c, (17) 



where b and are fermion annihilation and creation operators, respectively. 
They are defined ^Z5\ [2Sj by 



where we have chosen the convention a z \l) = — 11) and o~ z \0) = |0). 

Observe that H s is the XY-Hamiltonian which is well-known in solid 
states physics [25] and generates a quasi-free evolution of fermions. Since 
the so-called Jordan Wigner transformation, describing the fermion interpre- 
tation formally, is standard [27], we rephrase it only briefly. The subspace 
H, m can be reinterpreted as the space of m fermions moving without interac- 
tions in a 2m- dimensional state space. Then 7i m is simply the antisymmetric 
tensor product 



i=i 



b. ■= af^ 1 ® |0)(1| ® l® 2m -J, 
and satisfy the canonical anti-commutation relation 



b\bj + b 3 b\ = 6ijl, 




fit 
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The correspondence between these two pictures can be described easily: 
Given a basis state in the qubit register by a binary word with symbols 
1 at the positions j\, ■ ■ ■ ,j m - This state corresponds in the fermionic picture 
to 

(b'i)® \h)®---®\jm})~ 

where denotes here the jth canonical basis vector in C 2m . 

The restriction of H s to 7i m is in the fermionic picture given by 

H s = {S + 

+ 1® (S + &)®1---®1 
+ ... 

+ 1® 1® (S + S^), 
where S is the linear (non-unitary) shift acting on the basis states of C 2m via 

S\j) ■= \j + 1) 

for j < 2m and S\2m) = 0. This is because the term Yljbjbj+i shifts 
the fermion by one site. In other words, the time evolution in each tensor 
component in eq. (fTHj) is generated by the Hamiltonian S + 5', i.e., the 
adjacency matrix of the linear chain. This implies that the time evolution 
transfers annihilators and creators to linear combinations of annihilators and 
creators, respectively, when considered in the Heisenberg picture [27J. Let 

U t := ex P (-i(S + S^)t) (18) 

be the time evolution of one particle on a quantum walk on a chain and Uji- t 
its entries. Then the time evolution of the creation operator on site j is 

b] ~ E u ^ h \ » ( 19 ) 

l<2m 

i.e., it evolves into an operator creating a fermion which is in a superposition 
of different sites. For the annihilation operator at site j we obtain 

fri-X)^'* 6 '- (20) 

K2m 

To analyze the time evolution of relevant observables we will only make use 
of this formulation of the dynamics. 



20 



Time required for passing the circuit region 

The running time of our computer is given by the length of the time one has 
to wait until all atoms will have passed the interaction region R, i.e., a square 
of length k < n, almost with certainty. The dynamical evolution exp(—iH s t) 
is clearly quasiperiodic because the Hilbert space 7i m is finite dimensional. 
Therefore the atoms will always return to R. However, the essential idea of 
the ergodic quantum computer is that the probability of finding it outside of 
R is high if one measures at a random time instant. Later, we will therefore 
also consider the time average, before we estimate the time needed by the 
atoms to pass R for the first time. We found: 

Theorem 2 Let the circuit region be a square of length k. Then there is a 
time t G 0{k) such that the probability of finding all atoms outside the circuit 
region R is at least 1 — 12/ k given that the size n 3> k of the whole lattice is 
sufficiently large. 

The proof can be found in the appendix. Note that the probability of 
finding some atom in R could even be made smaller if the factor 4 in eq. ([2"oD 
in the proof would be replaced with some larger number. We conclude that 
the implementation time for our imprinted circuit is in O(k). 

Computing the time average state 

Now we want to show that also at a random time instant the probability of 
finding all atoms in the output region (i.e., the probability of finding at least 
k fermions on the left half) is close to 1. We found: 

Theorem 3 Let the circuit region be a square of length k = m/4 = (n— 1)/4. 
Then the probability of finding all atoms outside the circuit region tends to 1 
for n — > oo . 

This theorem is also proven in the appendix. One should maybe mention 
the following limitation of our model: in a lattice with finite length n the 
adiabatic approximation underlying the holonomic implementation is only 
true up to some error. Therefore the implemented transformation, given 
that the atom chain has passed the circuit region, does not exactly coincide 
with the desired quantum circuit. Instead, it depends to some extent on the 
time needed to pass the region. In the limit of averaging the time evolu- 
tion over an infinite time interval this error will increase. Roughly speaking 
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(and expressed in a too classical language), the atom chain has then traveled 
back and forth many times and the error can increase since passing the cir- 
cuit region backward would not necessarily invert the implemented circuit. 
However, it is a matter of the time scale on which the time average is taken 
whether this error is relevant, i.e., the degree to which the adiabatic approx- 
imation is justified determines the time scale on which the time average is 
described by our analysis. 

6 Remarks on the realization 

To judge whether it is likely that systems with Hamiltonians as above could 
be found in real systems would go beyond the scope of this article. Maybe 
one should rather ask which modifications are possible for our models that 
would make them more feasible. So far, the required spin-spin interactions 
are quite specific. 

To show that the required diagonal hopping is not a priori unphysical one 
could think of electrons on quantum dot arrays arranged as in Fig. The 
hopping along the rows does not require tunneling between distant dots even 
though it is the diagonal direction with respect to the square lattice. Spin- 
spin interactions would then only be needed in the direction of the dashed 
square lattice. 

Attractive interactions between particles on adjacent dots could, for in- 
stance, be achieved if the particles alternate between electrons and holes 
from row to row. Certainly, it will be difficult to find a system that combines 
the attractive interaction and the hopping with the spatially inhomogeneous 
spin-spin interactions. 

7 Imprinting a cellular automaton 

We have shown that nearest-neighbor interactions among qutrits located on 
a 2D square lattice can be designed in such a way that their autonomous time 
evolution simulates any desired quantum circuit. The Hamiltonian given here 
is significantly simpler than that one given in our previous paper [13] since 
the latter contained 10-qubit interactions. However, the advantage of the 
Hamiltonian in [13J is that it is programmable. The Hamiltonian given here 
contains the gates to be implemented already as hardware. In order to make 
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Figure 6: Schematic drawing of a dot array. The black circles indicate the particles 
carrying the qubits. Tunneling along the horizontal lines is possible since the 
distances are small. With respect to the dashed square lattice, the motion is along 
the diagonal direction. 
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Figure 7: Interaction regions that implement a cellular automaton according to a 
Margolus partition scheme. A pair of columns simulates one time step of the CA. 

a programmable quantum computer one could proceed as follows. Consider 
a universal quantum cellular automaton in one dimension [19] consisting of 
cells Cj with j G Z. Each cell contains a quantum system with Hilbert space 
C d . Using the so-called Margolus partition scheme [20j [191 HH1 [IS]> one time 
step of the CA consists of two parts. The first one is given by identical 
copies of a unitary U acting on all pairs (C 2 j_i,C 2 j), the second one by 
identical copies^ of V acting on all pairs (C^-, C^+i). In our scheme, we may 
represent each cell by the subspace of an appropriate number of adjacent 
logical qubits. Then we cover some part of our 2D lattice with a pattern of 
interaction regions that implement unitaries U and V as shown in Fig. [7J 
The number of pairs of columns in this pattern determines the number of 
time steps of the simulated CA. The translation symmetry of such a crystal 
would then, however, be described by huge unit cells. Nevertheless, our 
construction shows that translation invariant finite-range interactions among 
finite dimensional quantum systems can in principle implement a universally 
programmable quantum computer. 

2 To avoid technical problems with infinite tensor products of the form • • • <g> V <B> V <g> 
• • • one could use a C* -algebraic description and work with automorphisms on so-called 
quasilocal algebras [5D] . 



24 



8 Conclusions 



We have shown that nearest-neighbor interactions on 2-dimensional lattices 
can induce dynamical evolutions that are powerful enough to represent a 
programmable quantum computer. Even though we have not described a 
realistic implementation scheme we have given interactions which are not 
a priori unphysical. The main reason why it may be difficult to construct 
systems having exactly the Hamiltonians considered here is that the spatial 
homogeneity of the interactions follow a rather sophisticated law. 

To find even simpler Hamiltonians which can perform quantum comput- 
ing in a closed physical system is an interesting challenge for further research. 
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A Alternative model avoiding diagonal inter- 
actions 

The clocking scheme of our autonomous computer presented so far is some- 
what sophisticated in the sense that the hopping terms connect lattice sites 
lying diagonal with respect to the grid given by the possible atom positions. 
In other words, the rows defining the direction of atom propagation were di- 
agonal with respect to the square lattice structure. It would be more natural 
to have tunneling along rows or columns of the square lattice itself. 

We will therefore describe an alternative model for the propagation of an 
atom chain which has, however, the disadvantage, that we can only conjecture 
that the atoms move in forward direction with sufficient speed and will appear 
outside the interaction region with reasonable probability. We can, however, 
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Figure 8: Left: Possible steps of atoms when they start in the first column. Right: 
Possible atom configuration and possible motions 

prove, that the corresponding classical random walk would really implement 
the computation. 

The whole lattice consists of 2n x m sites. For reasons that will become 
clear later, we will color the sites such that we obtain a pattern similar to a 
chess board. As indicated in Fig. [8] the 2n atoms start in the first column. 
In contrast to our first model, the white sites are used, too. 

The rules for the propagation of atoms along the rows are as follows. An 
atom on a white site is allowed to move forward or backwards by one column 
if the atoms in the two adjacent rows are currently in the same column. An 
atom on a black site is only allowed to move one column forward or backwards 
if it is finally located in between two adjacent atoms in the same column. 
In the first column, only the atoms on the white sites are allowed to move 
forward. A possible configuration obtained by applying these rules is shown 
on the left of Fig. E 

By induction, we will argue that the configuration of atoms has always 
the following properties: 

1. Atoms located in adjacent rows are either in the same column or in 
adjacent columns. 

2. For every atom located on a white site there are always atoms on the 
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two adjacent black sites in the same column. The only exception are 
particles on white sites in the upper-most and lower-most rows which 
have certainly only one black neighbor. 

It is clear that these conditions are satisfied by the initial state. In order 
to prove that they are preserved we observe that an atom that is located on 
a black site can only move to a white site if the latter is adjacent to two 
occupied black sites. But then it will satisfy condition (2). When an atom 
on a white site moves to a black site, condition (1) is certainly preserved 
since (2) was true for its initial position. The atoms being on black sites on 
the boundary can only move to column j when the atom in the row below is 
already in column j. 

Now we present a Hamiltonian that induces these propagation rules. We 
start with 

2n to— 1 

K := ^2 °M 4,^+1 + h - c - > 
k=l 1=1 

which is a usual hopping term (as it appears in Hubbard models [28]) an- 
nihilating the atom at a certain position and creating it at the right or left 
neighboring site in the same row. Physically, one could think of an opti- 
cal lattice with periodic potentials generated by standing waves which result 
from the superposition of counter-propagating waves . To allow tunneling 
in horizontal direction and to avoid it along the vertical axis one could use 
high potential walls between rows and low potential walls between columns. 
This can be achieved by a superposition of high amplitude waves in vertical 
direction with low amplitude waves in horizontal direction. 

Now we modify the lattice such that there are two types of minima in 
the potential. The lower potential corresponds to the black sites and the 
higher potential to the white sites. By superposition of two lattices with 
wave length A and 2A one could also generate a chess-board like alternating 
potential (this is called a "superlattice" in |30j). 

We describe the alternating potential by introducing an additional term 

Ht := -2E N W , 

(k,l)EB 

where B denotes the set of black sites. Furthermore we introduce an attrac- 
tive interaction between atoms in adjacent rows which is only active between 
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pairs of atoms that are in the same column. The strength of the attractive 
interaction has strength E, i.e., we add a term 

2n— 1 m 

k=i 1=1 

This ensures that all the allowed atom configurations have the same energy: 
when an atom sitting on a white site moves to a black one the attractive 
interaction is switched off, this compensates the energy gap between the 
black and the white site. The same is true for an atom on a black field that 
moves to a white field since the attractive interaction is then switched on. For 
the upper-most and the lower-most rows we need to decrease the energy gap 
between black and white sites because there is only one attractive interaction 
compensating it. We achieve this by adding a potential 

H 3 := E ( Y, N V + E N ^j) ■ 

(l,l)eB (2n,l)eB 

Given the Hamiltonian Hi + H 2 + H 3 , one checks easily that an allowed 
forward or backward motion of an atom leads to a configuration with the 
same energy. Moreover, one can see that every step of an atom that violates 
the propagation rules would lead to a state with different energy. We observe 
furthermore that forbidden motions would always lead to a configuration with 
higher energy and never to lower energy. This additional feature provides the 
scheme with some thermo dynamical stability. 

Now we add the small perturbation K (E is again thought to be large) 
to the potentials and obtain 

H := H 1 + H 2 + H 3 + K . 

In analogy to the discussion of our first model we obtain an effective Hamil- 
tonian that is given by 

H eff :=PHP, 

where P is projecting onto the subspace of states having the same energy as 
the initial atom configuration. 

The effective Hamiltonian modifies the hopping terms such that they are 
multiplied by an additional number operator checking the position of the 
atoms in adjacent rows. In other words, each term of the form a^ k dj : k+i as 
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Figure 9: Spin-spin interactions for a logical one-qubit gate. 

well as those of the form a^aj fc+1 are only active if the two black fields are 
occupied which are either the current two neighbors of the considered atom 
(if the latter is located on a white site) or the future two neighbors (if it 
moves to a white site). Explicitly, we get therefore 

N N 

a Cg) a) + a) ® a + h.c. . 




The operators N of these 4-local terms act on black sites in the same column. 

To implement one-qubit gates in this scheme, we imprint interactions as 
shown in Fig. [91 The interactions V\, . . . , V\ are chosen as in Section [3j A 
decisive difference to our first model is that the hopping of the atoms does 
not turn Vj on immediately after Vj-i was turned off. Instead, the atoms 
pass always a configuration where all interactions are turned off. However, 
this is irrelevant for the holonomic scheme because the statement of Lemma [2] 
remains certainly true if the evolution is interspersed by time intervals where 
G(t) is completely switched off. 

To implement two-qubit gates, we design the spin-spin interactions as 
depicted in Fig.[TDl Now we have described how interactions can be imprinted 
that would implement the desired gates given that the whole atom chain 
passes indeed the circuit region. 

However, as already stated, it seems to be difficult to derive explicit 
formulas for the dynamical evolution. We will therefore only consider the 
corresponding classical random walk and argue that one will have at least 
probability 1/2 to find all particles outside the circuit region in the time 
average. Here the circuit region is defined by the first k columns of the 
chess-board where k is chosen such that the region contains all interaction 
stripes. The computation result is then obtained by measuring the spin of 
the particles after they have left the circuit region. 
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Figure 10: Spin-spin interactions for a logical two-qubit gate. Note that Uj and 
V are always switched on and off simultaneously since an atom can only visit a 
white field if the adjacent black fields in the same column are already occupied. 

The classical random walk would describe the physical dynamics in the 
limit of strongly incoherent atom hopping. The incoherent model is, of course 
inconsistent with our intention to study computation in a closed physical 
system. The motivation to study the classical walk is that the result gives 
us some small hope that the coherent walk behaves also as it is required for 
computation even though it may clearly be very different from the classical 
random walk. 

Let C be the set of allowed atom configurations on the 2n x m-chess- 
board and G be the graph with nodes C. Two nodes C\,C2 G C are adjacent 
if C2 can be reached from c\ by an allowed step of an atom. The probability 
distribution on the set of possible atom configurations at some time instant 
t is described by a vector p(t) := (pi(t), . . . ,Pi(t)) having the I nodes of G as 
indices. Then a continuous classical random walk on G is described by [31 J 

-p(t) = Lp{t) , 

where L is the graph Laplacian. Its entries are Ly := —1 for i ^ j if 
is adjacent, and Ljj = dj where dj is the degree of node j, i.e., its number 
of neighbors. Since G is connected, Ao := is a non-degenerate eigenvalue 
of L [31] and there is hence a unique stationary distribution. The latter is 
clearly invariant if we replace each atom configuration with the configuration 
obtained by reflecting it at the vertical symmetry axis of the chess-board. The 
latter is defined by a line between column 2/m and 1 + 2/m if m, where we 
have assumed that m is even for simplicity. 

Therefore the probability of finding at least one atom of the chain in the 
right half is at least 1/2. Since the chain does not tear this implies that the 
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column index of every atom is at least m/2 — 2n. Provided that k < m/2 — 2n 
we find with probability at least 1/2 all atoms outside the circuit region. 

B Proof of Lemma QJ 

We have a dominating Hamiltonian H pot and a weak perturbation K. In 
order to use a perturbation theorem of [8] we first have to rephrase some 
notations used in [H [32] . 

Given a Hamiltonian H pot and a small perturbation K, both acting on 
some finite dimensional Hilbert space 7i. Let H pot have a spectral gap A such 
that no eigenvalues lie between A_ = A* — A/2 and A + = A* + A/2 for some 
A* G R. Let be the eigenspace of H pot corresponding to the eigenvalues 
below A* and denote its complement by TL+. Denote the corresponding spec- 
tral projections by P±. For an arbitrary operator X, write X± T for P±XP T 
and X + for P + XP + . We define the self-energy operator for real- valued 

z by 

£_(*) := # po4 _ + + - , 

where G + , called the unperturbed Greens function (resolvent) in the physics 
literature, is defined by 

'.= zl + — Hpot^ . 

Using the above notations, we rephrase the following theorem which can be 
found as Theorem 8 in [8]. 

Theorem 4 Set H = H pot + K with \\K\\ < A/2. Let there be an effec- 
tive Hamiltonian H e ff with spectral width w e ff. We assume that H e ff = 
P_H e ffP_. Let < e < A and furthermore 

1. there is an r GB such that 

w e ff + 2e < r < A* , 

2. for all z e C such that \z\ < r, ||E_(2;) — H e ff\\ < e. 
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Then the restriction H<\„ of H to the eigenspaces with eigenvalues smaller 
than A* satisfies 



(21) 




(r - w eff - e)(r - w eff - 2e) ' 



We set A* := E/2 and A := E. Then the condition ||iT|| < A/2 is true 
since the expression in eq. (CQ) is a sum of less than n 2 terms of norm 1 and 
we have hence \\K\\ < n 2 < E due to the assumption E > n 6 . For the same 
reason, the spectral width w e ff of H e ff is less than n 2 . The operator P_ 
projects onto TL C and condition H e ff = P_H e jfP_ required by the theorem 
is satisfied by definition (see eq. (Hj)). We will now derive a bound on 

||£_(s) - H eff \\ 

in order to define an appropriate e and a corresponding constant r. 
Due to K = P_KP_ = H e ff and H pot _ = we have 

E_(z) - # e// = K_ + G + (l - K ++ G + )~ 1 K + _ 

The eigenvalues of H pot+ are E 1 , 2E, 3E, . . . according to the number of in- 
active interactions. Therefore the norm of 

G + (z) = (zl + + H pot+ y l 

is at most 2/E for all z with \z\ <E/2. With 

(1 - K++G+)- 1 = J2(K++G+T , 



we have 

||(l-^ ++ G + )- 1 ||<^ 
Hence we get 



n>0 



n>0 



2iiArii\" i 



E ) \-2\\K\\E 



\L_{z)-H eff \\ < \\K\\ 2 \\G, 
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1-2\\K\\/E ' 



For E > A\\K\\ we obtain with \\G+\\ < 2/E 

A r)4 

\\^-{z)-H eff \\<\\Kf-<A-=:e. 

We may choose r := \fE in order to fulfill w e ff + 2e <C r <C A*. We may 
now use inequality f[2"TT) and obtain 



\H E/2 -H eff \\ < 3(n" + l 
+ 



n 



2 



E-n 2 
4n 4 



VE-n 2 - 1)(VE 



where we have inserted the following results and definitions from above: 
w e ff < n 2 , A* = E/2, A = E, r = \[E. We have inserted e = 4n A /E 
only in the numerator of the second fraction and replaced it with 1 (as an 
upper bound) at the other places. Using E — n 2 — 1 > E/2 for E > n 6 and 
n > 10 and \[E - n 2 - 2 > \[E/2 we have 



He/2 ~ H eff \\ < An\ — + — - < 9 



2n 2 16n 



i n 3 



E E ~ Je 



l E/2 ~ n-eff 

where we have used lGn^/E < n 6 /E < n 3 /y/E. Thus, we have obtained the 
desired upper bound on the norm distance between H E / 2 and H e ff. 



C Proof of Theorem [2 

Let 

| J) := |0...01...1) e Tim (22) 

be the initial configuration when denoted in the qubit picture. Recall the 
beginning of section [5J where we have described the correspondence between 
clock states and binary words (as indicated by the sequence of symbols 1 and 
at the left of Figured]). 

Then all atoms are outside of R if and only if at least k symbols 1 have 
traveled from the right half of the chain to the left half, i.e., at least k 
fermions are contained in the left half of the interval of length 2m. We have 
therefore to solve a mixing problem of a "discrete free fermion gas" where all 
particles start in the right interval. First we show that after the time 0(k) 
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it is likely that at least k symbols 1 can be found on the left side. We define 
the observable 

m 

where Pj := b^bj is the projector on the upper state of qubit j. The observ- 
able M counts the number of symbols 1 on the left side. In Figure [2] this 
corresponds to the number of symbols 1 above the row in the middle. We 
will estimate the probability that less than k symbols 1 have moved to the 
left by the Chebyshev inequality. It states that for any random variable X 
we have 

P(\X-E(X)\>e)<^f-, (23) 

where E(X) and V(X) denote the expectation value and the variance of X, 
respectively. In the sequel we will consider the probability distribution of M 
as the distribution of such a classical random variable. 

Let \I t ) be the time evolved state after time t. Then the expectation value 
of M after the time t is given by 

m m 

EtW := ^(MPjlIt) = $>*I 6 ;W • ( 24 ) 

3=1 3=1 

Using the dynamics (TT^|) and (|2"UI) we get 

2m 

(I t \P 3 \I t ) = (I t \b]b J \I t ) = J2^tnMMW (25) 

i=i 

2m 2m 

= u ji;tUji;t = 



u ji;t\ 



l=m+l l=m+l 



where we have used that {I\b\bi\I) is 1 or depending on whether there 
the state \I) contains the symbol 1 at this position. Each term \uji- t \ 2 is 
the probability for a particle starting at site / to be found at site j when 
measured after the time t in a single particle quantum walk on a linear chain 
of length 2m. Since the time evolution of such a walk has been discussed 
in detail in the literature [33], we will only describe the implications for our 
model. 
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Consider a particle starting at site / with 

m + l<l<m + 4k. (26) 

In the notation of eq. (1221) such a fermion corresponds to one of the 4k 
leftmost symbols 1. The assumption n ^> 4k and hence m ^> 4k ensures 
that the time interval considered in the sequel is sufficiently small compared 
to the time to reach the boundaries of the lattice. Then reflections at the 
boundaries can be neglected. We have to wait only the time O(k) in order to 
achieve that the width of the wave function of a particle starting at a definite 
position is much larger than 4k (see [33J ) but still smaller than the size n 
of the whole computer. Then the probability of finding it on the left half is 
larger than 1/3. Recalling that this holds true for each I satisfying eq. (T26l) 
and that the number of expected fermions (= symbols "1") in the left half 
is given by the sum in eq. ff2^|) . we may choose a time instant t G O(k) such 
that the expectation value satisfies exactly 

E t (AT) = ^k. 

In order to estimate the variance 

V t (N) = E t (N 2 ) - (E t (N))\ (27) 
we observe that eq. (124"!) implies 

mu)f= (it\Pi\Wt\Pj\it) . (28) 

i ,j<m 

Moreover, we have 

E^N" 2 ) = ^2 (h\PiPj\h) = (It\PiPj\It) + E t (Af) , (29) 

where we have used 

m m 

Y,{It\PiPi\h) = Ys(UPi\h) = E t {N) , 

i=l i=l 

due to eq. (1241) . We rewrite one summand of the first term on the right of 
eq. ([29]) as 

2m 

{it\PiPi\h) = £ (I\b\b s b%\I) . (30) 

l,s,r,p=l 
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m<l<2m,s<m 

Ujs;t Ui s -t 

K/<2m, s<m l<l<m l<s<m 

|2 



The inner product can only be nonzero if annihilators meet creators, i.e., if 
either / = s and r = p or I = p and s = r or if all indices coincide. In the 
first case (including the third) the term is only non-vanishing for I = s > m 
and r = p > m since b\bi is the projection |1)(1| on qubit /. In the second 
case we must have I = p > m and s = r < m since b s b\ is the projection 
|0)(0| on qubit s. Hence eq. ( 1311 becomes 

^ Uu-t Uu-t U js - t Uj S ;t+ ^ Uu-tUjmUjs-tUis-t- (31) 

m<l, s<2m m<l<2m, s<m 

The first term coincides with (I t \Pi\It) (h\Pj\h) by eq. ( l25|) . We rewrite the 
second term as 

Uu-t Uji-t Ujs-t Uis-t 

<2m,s<rn 

^ Uu-t Ujl-t Ujs-t U is - ^2 U iht Ujl-t 22 
C2m, s<m 

= ~\ ^2 u H;tUji;t\ 

l<l<m 

where the last equality is due to 

uu-tUji-t = Sij , 

l<l<2m 

because U t is unitary (see eq. [IH|1 . Hence we have found 

(It\PiPj\It) < (I t \P t \It)(It\P 3 \It) Vi^j- (32) 
Combining (T2"T|) and fl32|) with fT28|) straightforward computation shows 

V t (AT) < (EtW - J2((W\It)) 2 ) < E t {J\f) . (33) 

Recall that we have chosen the time instant such that E t {M) = 4k/ 3 and 
hence V t {M) < 4k /3 by ineq. (1331) . Assume we would find less than k symbols 
1 on the left half. Then the random variable defined by Admeasurements 
would deviate at least k/3 from its expectation value. Hence we can apply 
eq. (T23D with V t {M) < 4k/ 3 and e 2 = k 2 /9 and find that the probability of 
finding less than k symbols in the left half is at most 12 /k. This completes 
the proof. 



36 



D Proof of Theorem 



3 



In analogy to the proof of Theorem [2] we will compute the expectation value 
and the variance of Af in the time average state and then use the Chebyshev 
inequality (1231) . 

The time average expectation value of the fermion number on the left 
half is given by averaging eq. (|24p over all t: 

E(M):=\im^[ E t (Af)dt = jjT \ f dt ) > ( 34 ) 

J ° l=m+l j=l ^° 

where the last equality follows from eq. (1251) . Recall that we interpret each 
summand for I — m + 1, . . . , 2m as the probability of finding a particle (that 
has started at site I) in the left half of the chain when measured at a random 
time instant. We will show that it is close to 1/2 up to an error in 0(1/ y/m). 
It is known [34] that the Hamiltonian S + S^ on C 2m generating the walk (i.e., 
the adjacency matrix of the "path graph P2m") has 2m different eigenvalues 

rir 

X r = 2 cos r = 1, . . . , 2m . 

2m +1 ' ' 

The eigenspaces of S + , are therefore one-dimensional and the rth eigen- 
vector is given by [35] 

2m 

\e r ) := v^^cos((j--)(r-l)^-)|j), 

3=1 

where the normalization factor is c = l/(2m) for r = 1 and c = 1/m for 
r ^ 1. 

Using these eigenvectors, we may write the time average density matrix 
of a particle that has started at position I as 

2m 

^2 \e r ) (e r \l) (l\e r ) (e T \ . 

r=l 

By evaluating the probability to be at position j using this state we may 
compute the time average in eq. (134"|) and obtain 

lim I [ T M 2 dt = f^\(j\e r )f\(e r \l)\\ (35) 
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Each eigenvector defines a probability distribution p r on {1, . . . , 2m} by 

PrU) ■= I0>r)| 2 = \ [l + cos ((j - ~)(r - l)£)] . (36) 

The spatial probability oscillations are given by waves with frequencies u r : = 
(r — \)ti /m. We will refer to those frequencies v that satisfy 

1 1 . 

2vr^= < i/ < 2tt(1 - 



m y/m 



as frequencies (they are neither close to nor close to 2n) and to the 
others as low frequencies. For a high frequency wave, the probability distri- 
bution j i— > p r (j) leads almost to equal probabilities for both halfs. This is 
seen from 

m 1 1 m 

X>r (j) = 2 + - E cos ^ ~ y) • 

j=i i=i 
The sum over the oscillating terms can be bounded from above by observing 

i E - \)\ ± i E eilu \ ^ G 0(v ^ } ' 

l<l<m l<l<m 1 1 

where we have used that v differs from both and 2n by fl(l/^ s /m) and the 
last inequality follows directly from the geometric sum formula 



71 1 _l_ a n + l 

^ H l-q 



Using eq. fl35l) and eq. fl36l) , the term in the bracket in eq. (JM]1 can be written 



as 

m 2m 



EE^O')l(er|0| 2 - (37) 

j=l r =l 

Neglecting the low frequency terms p r in this sum can only cause an error 
in 0(1/ y/m) due to |(e. r |j)| 2 G 0(l/m) taking into account that there are 
0(\fm) such terms. Hence the probability of finding a particle that has 
started at any site I in the left half is close to 1/2 up to an error in 0(1/ y/m). 
By summation over all / = m + 1, . . . , 2m we find E(M) — m/2 G 0(mj y/m). 
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To derive bounds on the variance of M in the time average state we use 
eq. (1291 and observe 

E(Af 2 ) = average((/ t |P i P J |/ t )) + E(M) (38) 

< Yl &vn^((h\Pi\It)(It\Pj\It))+E(Af), (39) 

where the last inequality is due to ineq. (1321) . We may rewrite the term in 
the big bracket in equation (1391) using eq. fl25l) into 

2m 2m 

(itlp^iitl^it) = M 2 Yl m 2 - ( 4 °) 

Z=m+1 /=m+l 

Recall that the coefficients uu- t are the matrix entries of a unitary that de- 
scribes a random walk on a linear chain. Note furthermore that the corre- 
sponding unitary group Ut is generated by a real-symmetric Hamiltonian and 
we have therefore uu- t = Uu :t . Thus, we are allowed to interpret \uu :t \ 2 as the 
probability for finding a particle at position I that has started at position i, 
which interchanges the original roles of i and I. This will be convenient in the 
sequel. The product of the sums on the right hand side of eq. (1401) can either 
be interpreted as arising from two independent walks of two particles or the 
walk of one particle in two dimensions. Explicitly, we consider a dynamical 
evolution in C 2m <g> C 2m generated by the Hamiltonian 

(S + ® 1 + 1 <g> (S + &). (41) 

Since the Hamiltonian (T4T]) is the adjacency graph of the square lattice 

{l,...,2m} x {I,. ..,2m}, 

we can consider the right hand of eq. (HOI) as the probability of finding a 
particle starting at \i,j) in a quantum walk on the square lattice in the 
"target quadrant" {m, . . . , 2m} x {m, . . . , 2m}. We want to prove that we 
have sufficient mixing in order to obtain 1/4 for the time average of term (T401) 
up to an error in 0(1/ \/m). We proceed similarly as for the one-dimensional 
walk with the essential difference that the spectrum of the Hamiltonian (T4T1) 
is degenerate since \e r , e p ) := |e r )<E>|e p ) and \e p , e r ) have the same eigenvalues. 
We denote the projection onto their span by P rp . The rank of P rp is 2 for 
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r p and 1 for r = p. Using these projections, we may compute the time 
average of (HOI) for fixed i,j by 

2m 2m 2m 2m 

Yl Y ( l > S \ P r,p\hj)(hj\ P r,p\ l , s ) = Yl ^'P^' S ) ' ( 42 ) 

l,s=m+l r,p=l l,s=n+l r,p=l 

where we have defined the coefficients 

2m 

d r,P ■= Y K^I^pN'-?)! 2 = (hj\Pr,p\i,j) 
l,s=l 

and probability measures q r ^ p on the square lattice by 

q r ,p(l,s) := -^—\(l,s\P r Ji,j)\ 2 . 

Defining the frequencies v r := (r — l)ir/n and v p := (p — ljn/n, the state 
vector given by normalizing P r , p \i, j) for a given r, p is a superposition of co- 
sine wave functions with frequency vector (u r , v p ) with another wave having 
the vector (z/ p , v r ). The corresponding probability distributions q r ^ contains 
then terms with frequencies 2u r , 2u p , v r — u p , v r + v p . We refer to a term as 
low frequency term whenever at least one of these frequencies is low. For 
each "high frequency distribution" q r ^ p the probability for the target quad- 
rant is 1/4 up to an error in 0(1/ y/rn) . This is seen in straightforward 
analogy to the corresponding argument for the one-dimensional walk. More- 
over, the contribution of the low frequency terms in (1421) can be bounded 
from above by the number of low frequency terms (which is in 0(m^ s /m)) 
times the maximal coefficient di >s of each term (which is in 0(l/m 2 ) due to 
| (i, j\e r , e p ) | G 0{l/m)). We conclude that the total contribution of all low 
frequency terms to (l4"2"j) is in 0(l/y/m). Hence the total sum in eq. ( T42l) is 
1/4 up to an error in 0(l/y/m). We conclude 

2 

average((/i|P,|/ t )(/ t |P J |/i)) = ^- + 0{m/^i) . 

Using eq. (1391) we conclude 

E(Af 2 ) = m 2 /4 + E(N) + 0(m/y/m) 
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and hence 

E(M 2 ) = (E(M)) 2 + 0(m) . 

Hence the variance E(N 2 ) — (E(N)) 2 is in 0(m) and the probability that a 
measurement of M leads to a result with less than m/4 converges to zero with 
0(l/m) by the Chebyshev inequality. Therefore, the probability of finding 
less than m/4 fermions on the left half (i.e., the atom chain has left the 
circuit region) tends to zero, too. 
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